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ABSTRACT 


For various applications in fluid dynamics, one can assume 
that the total temperature is constant. Therefore, the 
energy equation can be replaced ty an algebraic relation. 

The resulting sat of equations in the inviscid case is analyzed 
in this paper. It is shown that the system is strictly 
hyperbolic and well-posed for the initial value problem. 
Boundary conditions are described such that the linearized 
system is well posed. The Hopscotch method is investigated 
and numerical results are presented. 


This paper was prepared as a result of work performed under 
NASA Grant NCR 17-102-001 while the authors were in residence 
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1. INTRODUCTION 


The Navier-Stokes equations in two space dimensions con- 
tain four differential equations: the momentum equations, the 

continuity equation and the energy equation. For certain 
app * cations, the energy equation can be substituted by the 
as tion that the total temperature is constant, without 
much loss of accuracy, see (8]. The resulting system for two 
space dimensions is then in non-dimensionalized form: 


u t + uu x + 
v t * uv x + 


VU. 


♦ ip . 

px 


TpRe 


~ 14 (pv J -2 (pv„) +3 (p (v +uj)l 


X X 


y x 


vv + — p = •= — =r— (4(pv ) -2(pu ) + 3(p(v +u ) ) ] 
y p*y 3pRe y y H x y H x y' x‘ 


(1.1) P t + pu x + up x + p v_^ + vp y = 0 


2 2 

i = T + u + v (total temperature constant) 
p = RpT 


C 1 T 


3/2 


T+C. 


R, C^ f C 2 and the Reynolds number Re are given constants. 

In order to obtain good numerical solutions to the 
initial-boundary-value problem for the system when the 
Reynolds number is large, we must require that the inviscid 
equations (Re °°) are well posed. In this paper an analysis 
of the linearized version of the system is presented. The 
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characteristic speeds are no longer the same as for the com- 
plete inviscid Navier-Stokes equations, but the system is well 
pcsed for the pure initial value problem. The analysis of the 
mixed inltial-boundary-value problem shows that great care must 
be taken to obtain a well posed problem. If at a subsonic in- 
flow boundary, v and one of the variables u and p are 
specified, the equations are not well posed near the transonic 
point; v and a combination of u, v, p corresponding to an 
ingoing characteristic must in such a case be specified. 

Several numerical experiments have been done for this sys- 
tem by Rudy et al . [8). One method used by them was the 

Hopscotch scheme, see [3], [4]. The application of this 
scheme to the viscous terms is simplified if the function 
values at the middle point in the approximation of u , v 

A A A A 

are taken at time level n in both sweeps. It is shown in 
section 4 that this simplification introduces a stability 
limit on At. However, for high Reynolds numbers, it is more 
dissipative than the original method. 

2. THE PURE INITIAL-VALUE PROBLEM 

In this section we will first show that the linearized 
inviscid equations are strictly hyperbolic. After elimination 
of the variables p and T, the linearized system can be 
written in the form 

w. + Aw + Bw = 0 , 
t x y 

where 

a A <\ 

w - (U,V,F) 


9 
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A « 


(1-2R) u -2rv 
0 u 

P 0 


2 

c /P 


0 

u 


B 


V 

- 2 Rv 

0 


( 1-2R) v c /P 
P v 


2 2 2 
c* = R(l-u -V ), 

u,v,p now considered to be known functions. To prove strict 

hyperbol.’ city we must show chat the eigenvalues of Aw^ + Aw 2 

2 2 

arc real and distinct for all real , u> 2 with + w 2 <* 1. 

An easy calculation shows that these eigenvalues are given by 


X i = UO)^ + vu ) 2 


/ — 2 2 — 2 — 2 — 2 

X 2 3 = (l-R)(uii»j + vu 2 )f / R (uw 1 +vu) 2 ) +c (u^+u^) • 

All eigenvalues are obviously real, and since R ^ 0 , c ^ 0 , 

2 2 

they are also distinct for + w 2 = 1. 

The complete inviscid Navier-Stokes equations are 

symmetric hyperbolic, i.e., the corresponding matrices A and 

B can be symmetrized by the same similarity transformation, 

see [ 9 J . This is not the case for the system considered here. 

rr~2 — 7 

With the notation, a + = Ru±/R u +c , the eigenvectors of A 
arc the column vectors of 
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a + a _ 0 

S - 0 0 c 2 /R 

_ “P -P 2vp _ 

If A and B can be symmetrized by the same similarity 
tr« nsformation, then we can diagonalize A by an orthogonal 
transformation, and therefore B stays symmetric. 

Therefore S _1 AC is diagonal, and 



which is not symmetric. The eigenvectors in S can be 

permuted, but that does not effect the symmetry of S ^BS. 

Furthermore, each eigenvector can be scaled by different 

factors. This corresponds to a similarity transformation 

of S -1 BS by a diagonal matrix. It is easily shown that it 

is impossible to symmetrize S *BS by such a transformation, 

and therefore A and B cannot be symmetrized simultaneously. 

2 

The variable c defined by c = RT would for the non- 
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dimensionalizcd system correspond to the local speed of sound. 
However, the eigenvalues of A are u, (1-R) utc/l + Ru' /T. These 
are not the same as the eigenvalues of the corresponding A 
for the complete Navier-Stokes equations, since the latter 
are u, u i c. We note that these are obtained if terms of 
order R and smaller are neglected. (R is 1/7 for air.) 
Looking at the velocity component in the x-direction, a super- 
sonic state is most naturally defined as a state for which all 
the eigenvalues of A have the same sign. Therefore, we will 
define a state as subsonic if the condition 

(2.1a) ( 1-R) | u | < A 2 u 2 +c 2 

is fulfilled, and supersonic if the condition 

(2.1b) (1-R) | u | > /r 2 u 2 +c 2 

2 2 

is fulfilled. (2.1a) can also be written (l-2R)u < c , or 

equivalently 

(2.2) ^ u 2 +v 2 < 1. 

3. WELL POSEDNESS OF THE MIXED INITIAL- BOUNDARY VALUE PROBLEM 

In this section we will investigate the effect of the 
boundary conditions, and we assume the problem to be defined 
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on the domain 0^x< or >, - a »<y< 00 , O^t. We begin with 
a br'ef discussion of the one-dimensional case. 

Assume that the system w t + Aw x ■ 0 is transformed to 
diagonal form 


♦ t + A< t» x = 0, A « diag , X 2# X 3 ) . 

Then it is well known that the problem is well posed if the 
boundary conditions can be written on the form 


<fr X (0,t) = L<J> 1 1 (0 , t ) v g(t) , 

where (fr 1 contains those variables ^ which correspond 
to positive X^, and (J* 11 contains the remaining ones, see e.g. 
[6]. Using the same notation as in Section 2, we have in 
our case X^ = u , A 2 ■ u-a_ , X^ = u - a + , and 





7 


For supersonic flow (defined by (2.1b)), all three variables 

AAA 

U, V, R must be specified if u > 0 (inflow), and no 
boundary condition should be given if u < 0 (outflow) . For 
subsonic inflow, the boundary conditions must be such that 
♦ can be expressed in terms of $ ^ and an inhomo- 

geneous term. If we want to specify two of the physical 

A 

variables, we see that V must be one of them. Either one 

A A 

of U and R can then be chosen as the other specified 

( 2 ) 

variable, since for both of them 4> occurs with a nonzero 

A 

coefficient. For the same reason, either one of U and 

A 

P can be specified for subsonic outflow. 

The two-dimensional problem is much more difficult to 
analyze. The energy method does not work, since the system 
cannot be symmetrized, and therefore we must use the theory 
by Kreiss (5). Since part of the calculations are technically 
complicated, we give here a summary of the wellposedness 
proofs; Lho details are given in the Appendix. 

The wellposedness is determined by the behavior of the 
solutxons to the system 

A A 

(3.2) Aw v + (sI+iwB) w = 0, 

which is obtained by a Laplace transformation with respect 
to t and a Fourier transformation with respect to y. The 
solutions to (3.2) consist of components of the type e 1 i|». 

where is a solution of Det(C) = 0, C(tO^A*: + sI + iwB 

(or of the type x *e 1 ^ for multiple roots Kj). 
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With the notation a ■ s+iwv, the k^h are defined by: 

K l m ~ a/u 

(3.3) vk 2 j ♦ 2u( (l-Rja-RviwJiCj 3 * a 2 -2aR\ r ita)+c 2 w 2 “0, 

2 2 

where v « (l-2R)u -c . The corresponding vectors are 
given by 



uiui 


r- *■ 

*2 » 3 ^ a+,e 2, 3 U ^ 

*1 - 

a 

• *' 2.3 = 

iw(a +<2 3 U ) 


2 2 
2pR (u iw+va)/c 


p (u> 2 -k 2 , . { ) 


Note that a * ±uu> if and only if <2 = +a) and therefore 
<i * < 2 > In this case however, \p 2 is the nullvector, and we 
have only two linearly independent vectors ^ , i|^ . If 
<2 “ K 3 .it is also clear that there are only two linearly 
independent vectors ^ • 

A 

Let M be the class of vector functions w satisfying 

A 

(3.2) and with wf.L 2 (0, < > o ) for Re s > 0. The problem is well 

A 

posed if there is no nontrivial wcM satisfying the homo- 
geneous boundary conditions for any s with Rc s > 0. 
there is such a nontrivial solution for a purely imaginary s = s () 
then s o is called a generalized eigenvalue.) We will now inves 
tigate this condition for the different cases. 
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Supersonic Inflow 

In this case there are three linearly independent solutions 
belonging to I^tO,®), and we begin with 

Case It K \ ? k 2 ^ *3 ^ *1 

The general solution is 

w(x) * £e' # J ^ +ne K 2 V tj > 2 + 3^ iJ / 3 , 

where £#n#C are scalars. The condition for having nontrivial 

A 

solutions satisfying w( 0 ) = 0 is 
Det ^ 3 ) r 0 . 

A lengthy calculation (see the Appendix) shows that 
• (3.5) Det (4/ 1 ,»i» 2 ,^* 3 ) = v(k 2 -k 3 ) c 2 (u 2 ^ 2 -^ ) Z , 

9 

and it never vanishes because k 2 / and moreover, 

a = ± uo) only if ** < 2 • 

Case 2: < i ” K 2 ? *3 

The general solution is 

w(x) ■ +Cxij; 1 )e K l V ‘+n^ 3 e K 3 X , 

where 

C^l^l = - Af^ 


(3.6) 


\ 
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A 

The condition * 1**2 implies a * i uu and w(0) » 0 

implies When substituting this into (3.6), we obtain 

from the first and the last equation 

(3.7a) n(K 3 *w) 2 -P.U 

(3.7b) nliCjiw) 2 ■ l-pi-2Rpu (ui+v) /c^K , 

which contradict each other if £/0, n^O. 

Case 3: k i^ k 2 =ic 3 * 

The general solution is 

* . k ? x 

w(x> lX + (t 2 +nx^ 2 ) e # 

where 

(3.8) C(< 2 ) ^ 2 =-An «^ 2 • 

In the same way as above, we obtain two equations 
corresponding to (3.7) : 

(3.9a) iumpn=af; 


(3.9b) 


apn=-ui<i>£ . 
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2 2 2 

These lead to the case a «u w , which corresponds to 
k 1* ,,c 2' * nd t * 1 * 8 is 3 contradiction. 

Case 4; *c ^** c 2“* c 3 

In this case we must have a“±uw, k«+u>, and from (3.3) 
we get v<*-u [ ( 1-R) a-Rviw] . The imaginary part of the 
equation yields v*0 and tne real part yields u=0. 

This completes the proof of wellposedncss for supersonic 
inf lrv* 


Supersonic Outflow 

This case is trivial, since all k's have positive real 
part for Re s > 0, and no boundary conditions should be 
given . 


Subsonic Inflow 


In this case there are only two linearly independent 
solutions belonging to L^iO, 00 ). The general solution is 


(3.10) 


A ^ ^ A ^ 

w = £e Ij^+ne \p 2 , 


and 

P.ll) w = (J 1 +Cx^ 1 )e +UJX 


if 


if K^ = K 


2 
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where is defined by (3.6). 

A A 

With the assumption that U and V are specified, we get 
immediately the condition for a nontrivial solution 

(3.12) (uu) 2 +dK 2 ) (ci+k^’ 1 ) “ 0 if * 

Since by assumption cx+k^u^O, we obtain the equivalent con- 
2 

dition a*=-uu) With this a-value inserted into (3.3), we 

A 

obtain an equation for K^t^/iw: 

(3.13) (< 2 +l) (vk 2 -2Ruvk+u 2 ) = 0 . 

A A 

Since r = 4 i, a= 4 uu corresponds to k^=k 2 , the critical tc-values 
are 

(3.14) k * u (Rv±i4* 2 v 2 -v) /v 

The subsonic condition is crui/alent to v < 0; therefore, both 

A 

k-values are real. The corresponding a-values 
A ^ 

(3.15) a=aicj = - ■■ iw 

Rv±»^v 2 -v 

are purely imaginary. This means that the only nontrivial 

A 

solution corresponds to s Q = (a-v)iu>. With a=aiu)+6, 5>0, we 
want to see if the corresponding < satisfies Re<<0. In that 
case we have a generalized eigenvalue s=s Q , and the problem is 
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net well posed in the sense of Kreiss 1 5 ] . In the Appendix 
it is shown that Re k < 0 if 

(3.16) (1-R) (u 2 +v 2 ) > R. 

Unfortunately, this condition can be fulfilled even if the 
subsonic condition (2.1a) is satisfied, provided that 2R<1. 
Therefore, we have proved that the problem is not well posed 
for example with a transonic boundary, because there is 
always some part near the transonic point where (3.16) is 
fulfilled if v^O. 

It remains to treat the case with double roots 
when the solution h„.s the lorm (3.11). A straightforward 
calc lation (see the Appendix) shows that the condition for 
a nontrivial solution leads to the trivial case u=v=0. 

Another natural choice of boundary conditions would be 

A A 

the specif iciation of V and R. Assuming t ' Jie condition 

for a nontrivial solution is 

(3.17) c 2 a (<ij 2 -t< 2 )-2Riu)(u 2 iu>+VGi) (a+u< 2 ) = 0 • 

2 2 . 

Solving this equation for u> -tc 2 and inserting that into 
(3.3) gives (since a+uic 2 ^0) 


i 
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(3.18) (l-2R)uaK 2 «2Ru 2 uj 2 -a 2 . 

With the expression for defined by this equation inserted 
into (3.17), we obtain 

(3.19) (a 2 -u 2 <ii 2 ) (c 2 a 2 -4R 2 u 2 viu»(l-2R)a-4R 2 u 2 aj 2 (c 2 - ( 1-2R) u 2 ) ]^0. 


a=±uo) corresponds to te ^ *=tc 2 , so the critical a-values are 
given by the zeros of the second factor: 


(3.20) 


(• 


a=2Ruoj (Ruv ( 1-2R) it 


A 


2 __. 2„2 2 2 

vc - (1-2R) Ruv 


w 


From (3.18) it is easily seen that Rei< 2 > 0 if R**a>o 
and |a| <2Ru u>“. But t.iis inequality follows immediately 

from (3.19), where the magnitude of the constant term in the 
second factor equals c | ct j . He ice, the critical a-values 
(giving Ret^ 2 <0) must be imaginary, and this is the case when 


u 2 +(1+ g . 2 . ( - lyW 2 S ) V 2 , u 


A perturbation calculation shows that there is a generalized 
eigenvalue in a neighborhood of the transonic point. The 
above analysis show:; that wc must resort to the specification 
of the characteristic variables ^ and <t>* 2 ^ (see the 
definition (3.1).) Following the same lines as above, we 




i 


i 
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arrive for at t * ie equation (stated in terms of 

a this time) for a nontrivial solution: 

(3.21) (a 2 -u‘\> 2 ) (2Rviu>(u-a + )a + a-va 2 w 2 ) = 0 
2 

a =±ua) is ruled out by the assumption and the remaining 

critical a-value is, therefore, 

va + 

(3.22) a - - iu ‘ ai “ • 

A perturbation circulation shows that the condition Re k< 0 for 

A 

a=aiuj+6 , 6>0 , is equivalent to the supersonic condition. 

The multiple root case leads to no new restriction on 
u, v, and, therefore, we have a well posed problem. 

Subsonic Outflow 


For this case only one variable should be specified and 

a. KjX 

the solution to (3.2) always has the form = £e ij; 3 . Since 
a+< 3 u^0, the only possibility for a nontrivial solution with 

A 

U specified is < 3 = 0. This corresponds to the a-value 

/~2 2 2 A A 
a=(Fv±/c +R v )ia)=aiu). Substituting a=aiw+6 into (3.3), we 

obtain after dropping second order terms in 6 and < 
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K 




R ( Rv * /r v i- 


c*> 


6 


u<0 . 


Since R^v is always less than ( 1-R) /r^v^+c"^ if 2R < 1, we 
see that <>0 if 6>0, and wellposedness is proved. 

a 

For R specified , the condition for a nontrivial solution 
is <2 = 1(1). It is enough to investigate the case K 2 ~“ U) » uj>0 * 

From (3.3) we obtain the corresponding a-values <ij=uu),a,= 

( (1-2R) u+2Rvi ) a) which both have negative real parts. This 
proves the wellposedness. 

Let us finally mention that for a subsonic outflow 
boundary it is often difficult to specify accurate values for 
the variable u or p . In that case one could think of using 
numerical boundary conditions that approximate the condition 
3w r /3x r =0, r>0, i.e., vanishing derivatives of some order is 
assumed for all the variables. However, it is easily proved 
that this leads to a non wellposed problem for all r > 0. 

Since the derivative boundary condition applies also to the in- 
going characteristic variable, it is sufficient to study the 
scalar equation $ t +X<|> x =0, A>0. After transformation the 

A 

general solution is $=exp (-s/A ) , and the condition for a non- 

trivial solution becomes 
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(-s/A) r « 0. 


Therefore, s = 0 is a generalized eigenvalue, 


4 . THE HOPSCOTCH METHOD 


n 


Using the notation = u'jAx,nAt), the Hopscotch method 
(see (3), (4)) is defined by 


(4.1) u? + * = u 1 ? + Lu 1 ? , j+n even 

3 3 3 J 


(4.2) U j + 1 = U j + Lu j + ^» 3 +n oc *d 


where L is a difference operator. Assuming j+n even and 
combining equation (4.1) with 


u n - u 1 ? -1 + Lu n , 
3 3 3 


we obtain immediately the two equations 


(4.3) 


(4.4) 


n+1 « n n-1 

U 3 ' 2U 3 ' U 3 


u^ +1 = u^ 1 ♦ 2Lu ^ 


j+n even 


Similarly for j+n odd, we combine equation (4.2) with 


n n-1 . _ n-1 
u . = u . + Lu . 

3 3 3 
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to obtain 


..n+1 n-1 . _ . n+1 . n-1. 

■ Uj + L(Uj + ) . 


If L is defined by Lcr ■ (u " ^ L ~u } ) , the extra- 

polation formula (4.3) is valid with replaced by Lu^, and 


we get 


(4.5) u n+1 = u n ^ + 2Lu n , j+n odd. 

3 3 3 


Therefore, the hopscotjh scheme for the model equation u t »Au 
is equivalent to the Leap-frog scheme at every point, provided 
that the first time level for the latter one is generated Ly 
(4.1) and (4.2). The stability condition is the CFL-condition 


(4.6) || p (A) < 1 , 


where p (A) is the spectral radius of A. For the equation 


u =Au and with L defined by Lu n = A— 

XX 3 (Ax) 


, n _ n^ n . 

? <u j+l Jj-l 1 ' 


then, as pointed out in [3], the scheme is equivalent to the 
DuFort-Frankel scheme 


(4.7) u n+1 = u" -1 + A(u" ,-u n+1 -u n - l + u n ,) 

3 3 (ix) 2 3+1 3 3 3*1 
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with u* defined by (4.1), (4.2). Accordingly it is uncon- 

ditionally stable. In order to avoid the solving of a system 
of equations to obtain u n+ * at each time step, Au^ l+1 can be 
replaced by Au”. If only the steady state solution is wanted, 
the inconsistency with the time dependent problem thereby 
introduced if of no importance. However, we will show that a 
stability limit is imposed on At by this modification. 

For the scalar equation u^-ou , the modified method is 
defined by 


(4.8) u^ +1 =u ™+A (u^ +1 ~2u j +u j_^) , j+n even 

(4.9) u” +1 -u"+X(u^J-2u"+u"tJ) » j+n odd, 

where A *= — tq - j - . Proceeding along the same lines as above, 
(Ax j 

we derive the extrapolation formula corresponding to (4.3): 

(4.10) 2(1-A)u" = u" +1 + ( 1-2A ) u j _1 , j+n even. 

Using this we then obtain an equation containing function 

, n n-1 . 

values u j.j/ u j only: 

(4.11) u j + 1 =2A (1-A ) (u j + 1 +u j L ) + (1-X ) 2 u j -1 , j+n even. 

A trivial calculation shows that the von Neumann condition for 


(4.11) is 
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(4.12) 


X < 1, 


nt 1 


and the unconditional stability is lost. With u n 11# u n 

j±l j 

given, the equation (4.10) defines u^ assuming X < 1. 

We also want to investigate the dissipative properties 
for small X-values, which correspond to large Reynolds 
numbers in (1.1). The eigenvalues of the amplification matrix 
for the original Hopscotch method are denoted by z, and for 

A 

the modified version by z. They satisfy the equations 
(1+2X ) z 2 - 4Xaz - 1+2X - 0 


z 2 -4X (l-X)az- ( 1-2X ) 2 = 0, 


where a= cos(wAx). The solutions to these equations are, 
after expanding the square roots ano dropping 0(X S ) terms, 

z = 2Xa (1-X ) ± ( 1- 2X + 2X 2 ( 1+a 2 ) ) 

z = 2Xa (l-2X)+2X 2 oi± (1-2X + 2X 2 (l+a 2 )-2X 2 ) . 

A 

It is clear that if a^l, a^- 1 , then |z| < |zj. Accordingly, 
the dissipation is larger for the modified scheme, except 
for the lowest and the 1 ighest frequency, where both schemes 
have no damping at all. 
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There is also another way to look at this difference. 
From (4.8) and (4.9) we can for even j+n derive an equation 
on the form 


n+1 n-1. , n 

u. - Uj +2X(u jtl 




The last term represents the deviation from the original 
Hopscotch scheme, and it <s an approximation to 2AAtu^. 
Accordingly, the equation 


u t - T?T u xx 

is approximated by the modified scheme, and the dissipation 
coefficient is obviously larger for A < 1. 


5. NUMERICAL EXPERIMENTS 

The system (1.1) has been solved by Rudy et al. [8] 
using several numerical methods. For the high Reynolds 
numbers used there, we consider the system as a singular 
perturbation of the inviscid hyperbolic system. The choice 
of boundary conditions should therefore be made based on 
the analysis in section 3. 

When the experiments in (8] were made, only the one- 
dimensional analysis had been performed. Therefore, all the 
experiments were made with u and v specified at subsonic 


f I 1 I I 1 I 1 
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inflow boundaries. As we have seen in Section 3 , this ives 
rise to a non wellposed two-dimensional problem for transonic 
speeds at the boundary. However, for the actual boundary data, 
the condition ( 3 . 16 ) for a non wellposed problem was never 
satisfied at any grid point. In this section we will present 
results for another set of data, where condition ( 3 . 16 ) is ful- 
filled for the whole subsonic boundary. 

Figure 1 shows the computational domain and boundary data 
initially. In the neighborhood of the transonic point B the 
data were chosen in a way such that they were smooth on tno 
whole line AC. 

The boundary conditions on the line BC wore 


( 5 . 1 ) 


v nrl = £. 
oj J 


, e n n n+1 --..n n n+1 , , 2.n n+1 

( 5 . 2 ) -a .p . u . - 2 RV . p .v . + (c ) .p . =g. 

-O] o] 03 03*03 03 '03*03 J 3 


The outgoing characteristic variable t ^ was defined at the 
boundary by using linear extrapolation. At the upper 
boundary the analogous formulas were used. On AB every 
variable was specified, on EF linear extrapolation was used, 
and on the symmetry line A F we used the conditions 
Uy=py= 0 , v= 0 . This set of boundary conditions will be denoted 
by B . C . 1 . 

The scheme was also run with the subsonic inflow 

1 h : 


boundary condition ( 5 . 2 ) replaced by the condition u” * 1 
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and where was defined by extrapolation. This set will be 

denoted by B.C.2. In both cases the conditions 

|u n+1 -u n | < I0' 2 at|u n | 

(5.3) |v ntl -v"| < 10" 2 it|u n | 

I n+1 n i . , - ~ 2 . . | ni 
|p -p | <_ 10 At I p I 

were checked. For the Hopscotch fci enie an artificial viscosity 

2 2 

term approximating 0.1At((Ax) y +(Ay) Y ) was included in 

xx yy 

the equation for p. A 20x60 grid was used, and At=.019, 
Re=80,625. Figure 2 shows the pressure p after 500 time 
steps when conditions (5.3) were first fulfilled for B.C.l. 

It is seen that the pressure is far from being a constant in 
both c«..ses; B.C.2 produces very large and oscillating p-values. 
There are no oscillations for B.C.l. Figure 3 shows the 
pressure after 1300 steps. B.C.l gives a smooth and almost 
constant pressure profile. Figure 4 shows the results from 
B.C.2 after 1700 steps when the conditions (5.3) were first 
satisfied, and also after 2800 steps. Oscillations are still 
present in the whole subsonic region. 
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APPENDIX 
Proof of (3.4 ) 

The eigenvector is found directly by applying 
C(- » 0. 

However, since we don't have the explicit for^ of k 2 and 
the derivation of and ^3 b e m ore involved. Let 

i|>2 = ^ i » ^ 2 ' ^3 ^ 
and substitute it into 
C(te 2 )* 2 = 0 
to get with tc =< 2 

c 2 

[a+(l-2R)Ku]6 1 -2Rv<6 2 + ^-*63 * 0 

and 

c 2 

-2Ru.ia)6,+ (a+KU-2Rvia>) 6 ~ + — iui 6 , = 0 
1 2 p J 

Eliminating 6 2 from the last two equations yields 

2 

[ (a+(l-2R)+u) (a+KU-2Rviw) -4R uviu)K}<S^ 

_2 _2 

+ {— < (a+Ku-2Rviw) + — 2Rvtciu)}6, = 0. 

P p 3 

We use now (3.3) to get the expression 

2 

c 2 (k 2 -oj 2 ) 6 . = — <(a+<u) 6 , 

1 p 3 

and therefore, 

= ^(a+icu), = ~(tcu+a) , 5-j = - (< 2 -u> 2 ) . 
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Wc start with 


dot 


Derivation of Equation (3.5) 


— -(a+K-u) 
P 2 


— (a+< 7 u) 
P 3 


^•(a+< 2 u) 


•^(a+K^u) 
P 3 


2pR . 2 . . . 2 2 

-- ! ~ 2 (u iu+va) b) -K ^ 


2 2 
a ) -< 3 


Expanding the determinant around the first column yields 


2 2 2 2 

det (ip lf i|> 2 #i|> 3 ) = (< 2 -k 3 ) { (< 3 +< 2 ) [~2c co au+2Rau (viaw-u u )J 


(Al) + (< 3 -k 2 ) [-c 2 a> 2 u 2 -a 2 c 2 + 2Ru 2 (viaw-u 2 u 2 ) ) 


- c 2 u 2 <jj^-a 2 c 2 u ) 2 + 2Ra 2 (viau)-u 2 u> 2 ) } . 

We use now the characteristic equation (3.3) in order to 


(A2) k 2 +k 3 ■ ” 2u[ (l-R)a-2Rviu)]/v 


2 2 2 

(A3) = (a +c w -2Rviaia)/v . 


t 


r 
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After using (A2) and (A3) in (Al) we get 

det ( 4 ^, ^ 2 ^ 3 ) “ v(a 2 -w 2 u 2 ) 2 c 2 « 2 -< 3 ) 


Derivation of inequality (3. 16) 


With k and a defined by (3.14) and (3.15) respectively, 

A 

we insert a = aiw + 6 and k = Kiui + e into (3.3). Dropping 

2 2 

terms of order 6 , 6 e and e , we obtain 


e 


2 2 "2 

- 0 +c +R (u +u ) 



2 a 2 A 

where 0 = u il-R) + - Rva. Since u > 0, e is negative 

for 6 > 0 if and only if one of the inequalities 0 < 0 or 
(A4 ) 0 > c 2 +R (u 2 +a 2 ) 


2 

hold. With z = - v/(Rv) , the inequality 0 < 0 is equivalent 
to 


v ( + — ) > (l-R)u 2 , 0 < z < «. 

(1-/I+Z 2 1-/T+Z 

But since v < 0 and z > /1+z - 1, this inequality can never 
be fulfilled. 

The inequality (A4) is equivalent to 


(l-R)a 


Rva + v > 0 


f 
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which by the definition of a and z and after division by 
-v (>0) can be written in the form 


z(l-R) + 1 

(l±/ITzT 2 l±/I+z 


1 > 0 . 


2 

After multiplying by (lt/l+z) we obtain Rz + 1 < /1+z , 

which gives the inequality 


z < 


1-2R 

R 2 


Using the definition of z, v and c we get the inequality (3.16). 


The case k^=k 2 f° r subsonic inflow, U, V specified . 

Denoting the element.. of ^ by a^,a 2 ,a^, the last two 
equations of (3.6) are 

c 2 - 2 

(A5) -2Ruia^-2Rvia 2 + — ia^ = + u £ 

2Ru 2 

(A6) +a 1 +ia 2 = - — (uiivK 

For a^=a 2 =0, (A6) implies 5=0 for u^O, and from (A5) we then 

get a^=0, and there is no nontrivial solution. 


32 


A A 

The case with V and R specified 


A A 2 

With a"niu)*6 and K“*iu)+e, we get from (3.3) dropping 6 , 
2 

6 c and c order terms 


(A7) 


From 


A A 

u (1-R) K+a-Rv 

C — 7 

A 

-v|j - (1-R)a+Rv 
(3.20) we have 


6 

u * 


a « 2Ru(Ruv(l- 

A 

and the corresponding < 


A A 

( 1-2R) UKa = - 


2R) ±/(l-2R) 2 R 2 u 2 v 2 +vc 2 ) 
from (3.18) 

2 A 2 
2Ru^ - a . 


/c 


2 


Since |a| < Rv and sign k ■ - sign a, the numerator in (A7) 

A 

is negative if v > 0 and a > 0. Therefore e < 0 at the 
transonic point v = 0, and by continuity also in a neighbor- 
hood. 

T he case with Characteristic Variables Specified 

A A 

With a defined by (3.22), the corresponding k is given 
by 


(A8) 




va + . 2Ruv’ 

2 


2Rv (u-a + ) 


v 
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A A A 

With a"6+aico and K»c+tcioj we have (A7). Since sign ic sign 

A 

o *» - sign v it is sufficient to investigate the case v < 0. 
For wellposedness we want to prove e>0 which is equivalent to 


-vk/u - (l-R)a Rv > 0. 

A A 

Using the a - and < - values defined by (3.22) and (A7) we get 

2R 2 uv 2 (u-a ) 2 

(A9) (l-R)u|u-a + | < | v | + — a - -|--j . 

But it is easily seen that the subsonic condition is equivalent 
to the condition (l-R)u|u-a + | < |v|, and therefore (A9) is al- 
ways satisfied. 

For and w ^ th the notation $= (a 1 ,a 2 , a^JT , we get 

* ( i \ 

from the condition V=4- =0 on the boundary, that a^=0 and 

a l = ± 2 F y - (uiivK 

c 

C a 2 = ± (u pi + — ^ (uiiv)K • 

c 

( 2 ) 

The condition for a nontrivial solution satisfying $ =0 

on the boundary, becomes 


2Rv 



(-a +2Ru) 


2Ra u , n 2 2 

(1 =— + — r — i 


0 , 


and this cannot be satisfied since all terms of the imaginary 
part have the same sign. 



